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We point out that the transmission eigenvalue density and higher order 
correlation functions in chaotic cavities for an arbitrary number of incom- 
ing and outgoing leads (NijNz) are analytically known from the Jacobi 
ensemble of Random Matrix Theory. Using this result and a simple linear 
statistic, we give an exact and non-perturbative expression for moments of 
the form (A™) for m > — \N±— JV2 j — 1 and j3 = 2, thus improving the exist- 
ing results in the literature. Secondly, we offer an independent derivation 
of the average density and higher order correlation functions for f3 = 2, 4 
which does not make use of the orthogonal polynomials technique. This 
result may be relevant for an efficient numerical implementation avoiding 
determinants. 

PACS numbers: 73.23.-b,73.50.Td,05.45.Mt,73.63.Kv 



1. Introduction 

Conductance in mesoscopic systems is currently a very active area of 
research, both from the theoretical and the experimental point of view. In 
the scattering theory framework [1][2], for the case of a chaotic cavity with 
N\ and N2 channels in each of the two attached leads, the fluctuations of 
the transmission eigenvalues of the conductor are effectively provided by a 
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random matrix with appropriate symmetries [3] [4]. More specifically, the 
Dyson index (3 of the ensemble acquires the values 1 or 2 according to the 
presence or absence of time-reversal symmetry, or 4 in the case of spin-flip 
symmetry. 

Several quantities of interest for the experiments, such as the conduc- 
tance and the average shot noise, may be derived from the knowledge of the 
transmission eigenvalues {Aj}. Those are defined as the singular values of 
a transmission matrix t, which in turn is a N\ x N 2 off-diagonal block of 
a N x N unitary scattering matrix (where N = N% + N2) [5]. In the case 
of chaotic cavities considered below, {\j} are correlated random variables 
between and 1. 

Suppose that one is interested in computing the average shot noise (P), 
where: 

N 

P = P ^A P (1-A P ), N = min(N 1 ,N 2 ), (1) 
P =i 

Pq being a constant related to the physical properties of the conductor 
[5] [6] [7]. Until 2005, results for (P) were known only in the limiting cases 
JV ll2 > 1 [2] [8] [9] [10], Aq = N 2 = 1 [11] or few open channels [12]. Very 
recently, a compact form has been found for (P)/Po using two different 
methods, based on a semiclassical expansion [13] and on recurrence relations 
for the Selberg integral [14] (see also [15] for an alternative derivation). 
The latter nicely exploits the remarkably simple expression for the joint 
probability density (jpd) of transmission eigenvalues: 



N g 

p,(Ai, ...,\n)= a/-,- 1 n 1^ - n a ? < a,- < 1, 

j<k i=l 

(2) 

where the normalization constant is given by [14] [16]: 



N-i r (1 + 1 + jf) r U(\n 2 - m\ + 1) + if) r (1 + jf 
N p = J 



j=0 



f(|iV 2 -iVi| + l) + l + (iV + i-l)f 



r i + # r 



(3) 



A few comments about (2) are in order. The jpd in (2) is stated in [2] 
without proof and attributed to Brouwer. A formal proof has been given 
(using three different methods) by Forrester [17] in 2006, where the author 
also highlighted the connection with the jpd of the Jacobi ensemble of ran- 
dom matrices [16] [18]. In fact, one observes that the change of variables 
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yj = 1 — 2Xj brings (2) to the form: 

N 

Ppfa, ...,y N )= A/"" 1 J] ^-ykflKi-y^-^- 1 , -i < Vi < i, 

j<k i=l 

(4) 

allowing to use the machinery and results already known from Random 
Matrix Theory. 

In particular, the average density of transmission eigenvalues: 

N 

pp(X; N U N 2 ) = (V 5(X - A»)) =N d\ 2 ... dX N P {X, X 2 ...,X N ) 

(5) 

is of interest for computing linear statistics, i.e. observables of the form 
(tr/(ttt)>: 

(tr/(«t)>=/ dxpp(x;N u N 2 )f(x). (6) 
J o 

The moments of the form (A™) for a real number m can be computed in 
principle from the knowledge of the average density as: 

(Af) = [ dxx m pp(x;N 1 ,N 2 ), (7) 
Jo 

where the range for m is constrained by the convergence of the integral. 
The first two moments are directly related to the normalized conductance 
(G/Gq = (Ai)) thanks to the Landauer-Biittiker formula, and to the al- 
ready mentioned shot noise (P/Pq = (Xi) — {X\)). A refined semiclassical 
treatment of the former can be found in [19]. 

Surprisingly, the connection with the Jacobi ensemble has not been fully 
appreciated so far, with the consequence that the average spectral density 
Pb(X; N±, N 2 ) for finite and arbitrary number of open channels (Ni,N 2 ) is 
still deemed unknown (see e.g. [14] [20]). On the other hand, the density is 
known in the above mentioned limiting cases [2] [3] [4] [8] [12] [21]. 

In the mesoscopic literature the Jacobi ensemble is mentioned in the 
paper by Araujo and Macedo [12], where the authors derived the average 
density of transmission eigenvalues for a small number of open channels and 
(3 = 2 using an auxiliary non-linear sigma model. Their result reads: 

JV-l 

p 2 (X; N U N 2 < 11) = A" £ (2n + /i + 1){P^ (1 - 2A)} 2 (8) 

71=0 

where /u = f (|JV 2 -iVi| + 1) - 1 = |JV 2 -JVi|, N = min(JVi, JV 2 ) and P$? s) (y) 
is a Jacobi polynomial. 
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The authors state in [12]: 

...we believe (although we have no formal proof) that Eq. (8) 
is valid for arbitrary N\ and N2. This result is consistent with 
the random-matrix approach of Ref. [3] [4], which predicts for 
the same system a joint distribution of transmission eigenvalues 
given by the Jacobi ensemble, from which Eq. (8) can be derived. 
We have thus found independent evidence for the application of 
the Jacobi ensemble in this problem. 

However, the invoked references [3] [4] do not mention the Jacobi ensemble, 
and work out the only case Ni = N2. More precisely: 

1. Ref. [3] reports the jpd (2) restricted to the case N\ = N2 and (3 = 1, 2. 

For the case (3 = 2, the term A? ^ Vl ' +1 - ) 1 j n the jpd (2) then 
disappears, making the use of Legendre polynomials appropriate. For 
this subcase, the authors derive the average density and the 2-point 
function, and finally take the large N\ = N2 limit to get the smoothed 
macroscopic density /^(A; N\ = N2 S> 1) ~ N/iTy / A(l — A); 

2. Ref. [4] deals with all symmetry classes (3 = 1,2,4 and considers the 
two cases N\ = N2 S> 1 or N\ = N2 = 1. In the first subcase, the 
authors derive some quantities of interest with the use of a Coulomb 
gas approach after the change of variable Aj = 1/(1 + yi), yi € [0, 00). 

We wish to clarify that the average density of transmission eigenvalues 
for any N\ and N2 is exactly given by the density of the Jacobi ensemble, 
where the argument of the Jacobi polynomials is 1 — 2A (i.e. nothing but 
(8), for (3 = 2), and this result descends from the application of the standard 
Orthogonal Polynomial Technique [16] [22] to the (modified) jpd (4). In fact, 

the Jacobi polynomials P^ 1 ^ appearing in (8) are precisely the orthogonal 
polynomials over [—1,1] with respect to the weight (1 — y)^ in (4). The cases 
(3 = 1 and (3 = 4 are more complicated, but can be tackled in the same 
framework (see [23] and references therein). Also, n-th order correlation 
functions can be derived for all three symmetry classes [16] [23]. For example, 
for (3 = 2 one defines the kernel (see [16], Sections 5.7 and 19.1): 

N-l 

K N (x, y) = *^ V /2 E ( 2n + <" + l ) P n' Q) (! " 2x)P^'°) (1 - 2y), (9) 

n=0 

and the n-th order correlation function is written in terms of the (n x n) 
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determinant: 

Nl f 

p 2 (Ai, . . . , A n ) = ttt 77 / dX n+ ± . . . dAjvP/?(Ai, . . . , A at) 

= det[K N (Xj, Xk)]i<j,k<n (10) 

In particular, the average spectral density (1-point function) is exactly given 
by: 

N-l 

P2(X; iVi, N 2 ) = X^Y, ( 2n + V + 1){^ 0) (1 " 2A)} 2 (11) 

n=0 

extending the result (8) to an arbitrary number of open channels. In Ap- 
pendix B, we will show that for fi — > eq. (11) recovers the result by 
Baranger and Mello [3]. 

The purpose of this paper is thus twofold: 

• having clarified the role of the Jacobi ensemble, and the known results 
for its spectral density for arbitrary N± and N 2 , we give a closed form 
expression for moments of the form {Xf") for m > —\N 2 — N\\ — l and 
(5 = 2 through a simple integration over the average density (linear 
statistic). The formula is exact and non-perturbative, and extends 
previous results in the literature (Section 2). 

• exploiting a less known result by Kaneko, we give an alternative rep- 
resentation for the average density and higher order correlation func- 
tions for j3 = 2,4 in terms of hypergeometric functions of a matrix 
argument. Thanks to a recent algorithmic progress, this result may 
prove useful for a numerical implementation which avoids the use of 
determinants (or quaternion determinants) (Section 3). 

2. A closed form expression for moments 

For simplicity, we consider again the (3 = 2 case as in [12]. The moments 
(A™) can be computed as a simple linear statistic on the transmission eigen- 
values: 

(X?)= [ dxx m p 2 (x;N 1 ,N 2 ) (12) 
Jo 

Known results about (A] 71 ) include: 

1. Approximate evaluation for all positive integer m (but valid in the 
regime N±,N 2 ^> 1) [24] [25]; also, the generating function for all mo- 
ments in this limit was first computed in [21]. 



G 
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2. Exact evaluation (valid for all Ni,N2,f3) but only up to m = 4 (see 
[20] and references therein). 

Assuming N\ < N 2 without loss of generality, we can use (11) and (13) with 
Ni= N and N 2 = n + N: 

N-i 1 
(\f) = V (2n + n + 1) / dxx^ +m P^°\l - 2x)P^' 0) (l - 2s) (13) 

After the change of variables x = (1 — t)/2 and the definition of Jacobi 
polynomials as: 

P (M,0)m- 1 (-w)fc(^ + n + 1 )fc(/ i + fe + l )n~k f l-t \ k hA s 

^ fc! I 2 J (Mj 

fc=o v 7 

(where (a)^ = T(a + k)/T(a) is a Pochhammer symbol), we obtain: 

,, m > _ _J v-^ 1 2n + /x + 1 {-n) k {n + n + l) fc (/i + k + l) n _ fc 

^ 1 ' ~ 2^+™+! ^ n! ^ 2 fc fc! X 

n=0 fc=0 
X ^ (it (1 - t)^ +m+fc P^'°)(t) (15) 

The integral above can be computed for m > —fj, — 1 ([26], formula 7.391.2) 
in terms of a hypergeometric function 3i 7 2(— n, /i + n+l,/i + m + /c + l;/i + 
1, fj, + m + k + 2;l). Since the first argument is a negative integer, the series 
gets truncated to give eventually 1 : 



N-l 



{x ?) = y (2n + „ + 1) y — 9(kMi) (16) 



where: 



Despite lacking the aesthetic appeal of subcases already considered in 
the literature [24] [20], formula (16) is nevertheless valid for any (Ni,N2) 
and m > —fi — 1, and is fully non-perturbative. After implementing (16) in 
Mathematica®, one can check by direct inspection that: 



1 We are grateful to Marcel Novaes for suggesting significant simplifications in (16). 
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1. the formula (16) agrees with the approximate result in [24] (valid for 
N U N 2 » 1): 

w-<*+w±£:iy-^$^y (is) 

where c p = ^(J) (see Table 1). 

2. The shot noise power (P)/Pq, defined as (Ai) — (Aj), can be extracted 
from (16). Thanks to multiple cancellations, the result can be cast in 
the very simple form: 

(P) = N^ + Nf 

P (fi + 2N- l)(/i + 2iV)( M + 27V+ 1) 1 ; 

which agrees with the known exact result [13] [14] (see also eq. (30) 
below). 

3. The average conductance (G)/Gq = (Ai) from (16) can be brought to 
the simple form: 

(G) _ N(y + N) 

G ~ n + 2N [2U) 
which agrees with the known result [3]. 



A* 


N 


m 


Exact (16) 


Approximate (18) 


4 


57 


3 


18.4240 


18.4248 


4 


87 


7 


18.637 


18.638 


12 


47 


19 


6.7672 


6.77002 


15 


57 


29 


6.67909 


6.68199 


25 


75 


59 


6.34394 


6.34704 



Tabic 1. Comparison between the moments (A™) computed by Novaes [24] and 
our exact derivation (16). Note that the normalization J Q dx Pff(x; N\, N%) = Ni 
implies that the moments are not constrained between and 1. 



3. A second derivation of the average density and higher-order 

correlation functions 

In this section, we will derive an alternative expression for the average 
density of transmission eigenvalues and higher order correlation functions 
for finite N% and A2 and (3 = 2,4, starting from the jpd (2). Exploiting 
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a variant of the Selberg integral evaluated by Kaneko [27], all correlation 
functions can be expressed in terms of a hyper geometric function of a matrix 
argument, instead of a determinant of a kernel as in (10) (for (3 = 2). 
Consider the joint probability density of transmission eigenvalues (2): 

N 

p^ai, ...,x N )= n I** - m p n xj m - Nii+i) -\ o < a, < i, 

(21) 

where N = min(iVi, N2), ft = 1, 2, 4 and the normalization constant is given 
by (3). 

The density of eigenvalues is given by the following multiple integral: 



Pi 3 (X 1 ;N 1 ,N 2 ) =N f ... f d\ 2 ...d\ N P^Ai, ...,X N ), 
Jo Jo 

such that the normalization J" 1 pp(x; N±, N2)dx = N holds (where again 
N = mm(N 1 ,N 2 )). 

It turns out that the integral above can be evaluated without the use 
of the Orthogonal Polynomial technique, which would lead to the formula 
(11), if one resorts to the following extension of Selberg integral given by 
Kaneko [27]: 

„ n n 

/ u dx jH4^- x ^ 2 n (xi-tk)n\xj-xkf= 

J[0,l] n j =l j =l l<i<n j<k 

l<k<m 

Ci 2 pf /2) f-n, +£ 2 + m + 1) + n - 1; ~(4 + m); {h, . . . , t m }^j 

(23) 

where C\ is a known constant and 2^1 is a hyper geometric function of a 
matrix argument. Details about these objects are provided in the appendix. 
From (22), one has: 

IfljVa-jVxI+iJ-l 

pp(Xi;Ni,N 2 ) = 77 / dX 2 ...dX N Y\\X j -X k f 

N P V[ ,i]*-i f£ 

x J| A f(^-W)-\ (24 ) 

i=2 

Now, the Vandermonde coupling can be decomposed as: 

N 

J{\X 3 -X k f= n ^-X^Jl^-X/, (25) 

j<k j<k,j=2 j=2 
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and, for (3 = 2, 4 the absolute value in all products is immaterial. Hence: 



ppiXuNuN*) 



^fd^-MI+D-i 



/ <iA 2 . . .d\]\f | Xj — X 

./fO.lF- 1 , ,7: „ 



'[0,1] 

N 



j<k,j=2 



(26) 



i=2 



3=2 



Comparing (26) and (23), after the following identification: 



n 


= N-1 




h 


= f(|iV 2 -iVi 


+ 1)-1 




= 




tk 


= Ai VA; = 


1, . . . ,m 


m 


= 





one eventually obtains: 



pp(X 1 ;N 1 ,N 2 



iV Ci Af 



f (|JV 2 -JVi|+l)-l 



x 2 F| p/z; (1 - AT, \N 2 - Nx\ + N + 1; \N 2 - + 3 - 2/(3; X^) , 

(27) 

where we have introduced a customary matrix notation in the last argument 
of the hyper geometric function. Note that the result (27) is still formally 
valid for any even integer (3. 

We also observe that higher order correlation functions can be easily 
written down, exploiting the very same eq. (23). For example, the two- 

( 2 \ 

point function p\ (Ai, A 2 ; Ni, N 2 ) can be written (ignoring prefactors) as: 
pf(X 1 ,X2;N l ,N 2 ) oc (A 1 A 2 )Kl^-^l+ 1 )- 1 |A 2 - A^ 

x / axs-.-dx^xf^-^' 1 n \^-hf 

J l u > l l i=3 j<k,j=3 

N 

(28) 



N 



J=3 j=3 



and the (N — 2)-fold integral is again of the same form as (23) for the 
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following values of parameters: 



n 


= N -2 






h 


= lm 


-iVi 


+ 1)-1 


h 


= 






H: 


= Ax 


\/k = 


1,...,/? 


H: 


= A 2 


Vk = 


+ 1,.. 


m 


= 2(3 







Hence, this time the matrix argument of the hypergeometric function is 

X^ 2 ^ := diag Ai, . . . , Ai, A2, • • • , A2 I Note that the 2-point correlation 
\ p times /? times / 

(2) 

function p^ (X\, X2; N\, N2) is manifestly symmetric in the two arguments 
as it should, due to the symmetry of Jack polynomials (see appendix). It 
is worth mentioning that higher order correlation functions can be written 
down easily along the same lines. 

Thanks to a very efficient Matlab® implementation of this kind of 
hypergeometric functions by Koev and Edelman [28], the density itself, lin- 
ear statistics (one-dimensional integrals over the density) and n-th order 
correlation functions can be numerically tackled very easily. In particu- 
lar, these results entirely avoid the use of (quaternion) determinants and 
(skew-)orthogonal polynomials which would arise from the canonical RMT 
treatment and can get computationally demanding for high iVi2 and n. 
Conversely, the computational complexity of the algorithm in [28] is only 
linear in the size of the matrix argument {(3n). 

In the following, we shall provide some plots of the spectral density for 
different numbers of incoming and outgoing leads, and (5 = 2 (Fig. 1). The 
agreement between the two alternative formulas (11) and (27) is excellent. 

As a final check, we also numerically compute the prototype of linear 
statistics, i.e. the (normalized) average shot noise power (P)/Pq (see (1)), 
defined as: 



(P)/P = [ dxp /3 (x;N 1 ,N 2 )x(l-x) (29) 
Jo 

where pp{x\ N±, N2) is taken from (27). The result has to agree with the 
analytical expression [13] [14]: 

CP> = jVi(iVi-l + 2/ / 9)jV 2 (jV2-l + 2//3) 
i^o (N -2 + 2/ (3)(N -l + 2/(3)(N -1 + 4/(3) 

where N = Ni + N2- We compare in Table 2 the theoretical result (30) with 
the numerical integration of (29), obtained in MATLAB® with a standard 
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integration routine. The agreement we found is excellent, thus confirming 
the correctness of (27). 



35 



30 





N. 


= 6; N 2 = 


8 (Jacobi) 




N 


= 6; N 2 = 


8 (hypergeometric) 


o 


N 


= 2;N 2 = 


4 (Jacobi) 




N, 


= 2;N 2 = 


4 (hypergeometric) 


if 


N i 


= 10;N 2 


= 14 (Jacobi) 




N. 


= 10;N 2 


= 1 4 (hypergeometric) 



25 - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Fig. 1. Density of transmission eigenvalues for (3 = 2 and different values for the 
pair (iVijA^). The plot symbols are used for RMT formula (11), whereas solid 
lines represent the alternative formula (27). 



iVi 


N 2 


P 


Theory 


Numerical 


4 


7 


2 


0.5939393 


0.5939393 


8 


11 


2 


1.1321637 


1.1321639 


3 


9 


2 


0.4248251 


0.4248251 


4 


7 


4 


0.5805422 


0.5805424 


3 


5 


4 


0.4326923 


0.4326923 



Table 2. Comparison between the theoretical expression for the average shot-noise 
power (30) and the numerical integration of (29), for different values of N±, N2 and 
0. 
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4. Conclusions 

We have clarified that the average density of transmission eigenvalues 
for chaotic cavities is analytically known from the Jacobi ensemble of ran- 
dom matrices, as well as all higher-order correlation functions. The known 
formula for the average density coincides with the one derived with a map- 
ping to a non-linear sigma model by Araiijo and Macedo, but is rigorously 
valid for any number of open channels. With an elementary integration over 
this density, we derived a general formula for the moments (A™), which is 
non-perturbative and valid for arbitrary large m and (3 = 2. Also, a second 
derivation is offered for the spectral density and higher order correlation 
functions, which does not make use of orthogonal polynomials or deter- 
minantal structures. Thanks to a recent algorithmic progress, this result, 
exploiting a hyper geometric function of a matrix argument, may be numer- 
ically easier to implement than high-order (quaternion) determinants. All 
the results are consistent with numerical checks and known formulas in the 
literature. 
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Appendix A 

Hyper geometric function of a matrix argument 

Following Kaneko [27], we first report the definition of the constant C\ 
appearing in (23): 

Ci :=5 n , (4 + m^2,/3), (A.l) 

where: 

q , fl „ A . A r(»f + i)r( yi + i + (i-i)f)r( y2 + i + (i-i)f) 
s n ,o( yi ,y 2 ,z) -11 r( i + i)r (2/1 + y2 + 2 + (n + ,-2)|) " 

(A.2) 

The hypergeometric function of a matrix argument [29] takes a sym- 
metric matrix (m x m) X as input and provides a real number as output. 
It is defined as a series of Jack functions of parameter j3, which generalize 
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the Schur function, the zonal polynomial and the quaternion zonal polyno- 
mial to which they reduce for j3 = 1,2,4 respectively. Given a partition 
n of an integer k, i.e. a set of integers k% > n 2 > • • • > such that 
| | = «i + «2 + - . . = fc, and a matrix X, the Jack function d^(X) is 
a symmetric and homogeneous polynomial of degree \k\ in the eigenvalues 

X\t ■ ■ ■ > of X. 

The hypergeometric function is defined as: 

pF p( ai , ...,a p -M,...,b q ;X):=J2Y: SP^^^, 

(A.3) 

where the symbol k \- k means that k is a partition of k and (a)J^ = 
ri(j j)e K ( a ~~ ^ + J — is a generalized Pochhammer symbol. 

The series (A.3) converges for any X if p < q; it converges if maxj \xi\ < 1 
and p = q + 1; and diverges if p > q + 1, unless it terminates [28] [29] [18]. 

Appendix B 

77ie // — > /irm£ o/ £/ie spectral density. 

In the case N\ = N 2 = N and = 2, the average spectral density was 
computed exactly in [3] as: 

N 2 

p 2 (X;N 1 = N 2 = N) = — — — {P%(a) - 2aP N (a)P N . l {a) + P 2 N ^{a)} 

(B.l) 

where a = 2A — 1 and Pn( x ) is a Legendre polynomial. 

This case corresponds to fj, — > in eq. (11). In this appendix, we show 
explicitly how to get from (11) to (B.l). 

First, we remark that the identity between Jacobi and Legendre polyno- 
mials P^'°\x) = P n {x) holds [formula 22.5.35 in [30]]. Hence, in the case 
/i->0we have from (11): 

N-l 

P2 (X; N 1 = N 2 = N) = Y, (2« + l){Pn(l - 2A)} 2 (B.2) 

n=0 

Next, we use the Christoffel-Darboux formula for Legendre polynomials [see 
formula 22.12.1 in [30]] at equal arguments x = y = 1 — 2A: 

JV-l 

(2n + l)Pn(v) = N[P' N {y)P N ^{y) - iViCvWlO] (B.3) 

n=0 



14 



Transmission printed on March 10, 2008 



Then, we exploit the differential relation [22.8.5 in [30]]: 

(1 - y 2 )P^(y) = -nyP n (y) + nP n _i(y) (B.4) 

to get: 

p 2 (X;N 1 =N 2 = N) = -JL-jlNPfr^y) - yP N (y)P N ^(y) - (N - l)x 

x P N {y)P N _ 2 (y)] (B.5) 
y=l-2X 

Thanks to the recurrence relation [22.7.10 in [30]], we obtain the following 
identity for Pn-2(v)'- 

P N -2(y) = JfrjK 2N ~ 1)vPn-x(v) ~ NP N (y)} (B.6) 

which is then substituted into (B.5). Eventually, given that y = 1 — 2A and 
the Legendre polynomials have the same parity of their index, we obtain: 

/V 2 

p 2 (A; N 1 =N 2 = N) = — — — [P^(2A - 1) - 2(2A - 1)P N (2X - l)x 

xP JV _i(2A-l) + P^_ 1 (2A-l)] (B.7) 
in complete agreement with (B.l). 
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